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We study the inflated phase of two dimensional lattice polygons with fixed perimeter A'' and 
variable area, associating a weight exp[py4 — Jh] to a polygon with area A and h bends. For convex 



and column-convex polygons, we show that {A)/A^^ 



1 - K{J)/p' + Oip-"), where p = pN -> 



1, and p < 1. The constant K{J) is found to be the same for both types of polygons. We 
argue that self-avoiding polygons should exhibit the same asymptotic behavior. For self-avoiding 
polygons, our predictions are in good agreement with exact enumeration data for J = Q and Monte 
Carlo simulations for J 7^ 0. We also study polygons where self-intersections are allowed, verifying 
numerically that the asymptotic behavior described above continues to hold. 



PACS numbers: 05.50.-|-q,02.10.Ox,05.70.-a 



I. INTRODUCTION 

The study of lattice polygons weighted by area and 
perimeter is a central problem in lattice statistics and 
combinatorics. Lattice polygons have been used to model 
vesicles P, cell membranes 0, emulsions Q, poly- 
mers and percolation clusters 0. In several cases, 
exact generating functions for classes of such polygons 
have been obtained. A survey of different kinds of lattice 
polygons and a review of related results can be found in 
Refs. 0,i. 

In this paper, we study the asymptotic behaviour of 
the area enclosed by inflated polygons of fixed perime- 
ter. We calculate the area for two special cases of lattice 
polygons - convex and column-convex lattice polygons. 
We then conjecture the appropriate form for the area of 
self-avoiding polygons in the inflated phase. 

We first summarize known results for the problem of 
pressurized polygons, based on the generating function 



(1) 



A, TV 



where Cm{A) is the number of self-avoiding polygons of 
perimeter N and area A, weighted by a chemical poten- 
tial ^. Here p is the pressure which couples to the area A. 
Exact solutions exist for G{u,p) when Cn{A) is restricted 
to convex polygons [13, [HI or to column-convex poly- 
gons [13]. However, a general solution for self-avoiding 
polygons is unavailable. Exact enumeration results for 
self-avoiding polygons exist for all up to = 90 and 
for all A for these values of N [l^l . A transition at p = 
separates a branched polymer phase when p < (for /i 
sufficiently small) from an inflated phase when p > 0. At 
p = 0, the problem reduces to that of the enumeration of 
self-avoiding polygons. The scaling function describing 
the scaling behavior (for p < 0) near the tricritical point 
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p = and /.J = K^^, where n is the growth constant for 
self- avoiding polygons, is also known exactly [l3. [isl. [l6| . 

Less is known about the inflated phase obtained for 
positive pressures p > Q. However, this phase is of phys- 
ical interest in the case of two-dimensional vesicles, or 
equivalently pressurized ring polymers [ij, [ia,[ii,[ii,[a. 
In the calculations described in this paper, we consider 
the partition function 



pA-Jb 



P>0, 



(2) 



A.b 



where Cpf{A,b) is the number of self-avoiding polygons 
of area A with b bends. A bending energy cost J per 
bend is introduced to incorporate semi-flexibility. 

Some rigorous results exist for Zff{p,0) when p > 0. 
Ref. [m proved that 

Zjv(p,0) = A(p)ef^'/i6(l + 0(p^)) as A^ -> 00, (3) 

for some < p < I, with A{p) some function of p. This 
result holds in the limit where N 00 keeping p fixed. 
However, as we argue in Sec. [Hi the relevant scaling limit 
in the inflated regime is p ^ 0, A^ ^ 00 keeping p = pN 
finite. 

The central result of this paper is then the following: 
In this limit, we show that for both convex and column- 
convex polygons, the area is given by. 
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^ 327r2 64^. , 
1 - -tt;^ + -35-L12 (1 - a) 
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where Li2 is the dilogarithm function 



+ 0(e-P/«), (4) 



(5) 



and, a = . We argue that this result should also 
extend to the self-avoiding case and test this conjecture 
numerically. 

The paper is organized as follows. In Sec. [Til we present 
a justification of the scaling limit we consider using a sim- 
ply Flory-type argument. Sections. IIIII and IIVI contain 
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the calculation of the area for convex and column-convex 
polygons respectively. Section |V] contains the numeri- 
cal analysis of self-avoiding polygons and self intersecting 
polygons. A brief summary of our results and conclusions 
is presented in Sec. IVII 



II. 



SCALING IN THE INFLATED REGIME 



The inflated regime of self-intersecting pressurized 
polygons has been well studied in the continuum [TtI . [isl . 
I19I . [20|. In this case the appropriate scaling variable is 
obtained by scaling the thermodynamic pressure with the 
system size, taking p ^ 0, N ^ 00 keeping p ~ pN fi- 
nite. A typical configuration in the inflated phase has no 
self-intersections. Thus, we expect that the above scaling 
should also hold for self-avoiding polygons. 

To motivate our choice of the scaling variable, we 
shall follow the Flory-type scaling analysis developed in 
Ref. [23. The free energy in the inflated phase consists 
of two terms. These describe the contribution from the 
pressure differential and the stretching free energy of the 
closed self-avoiding walk. The pressure contribution is 



pressure 



= -pA « -pR^ 



(6) 



where it is assumed that the area scales as the square of 
the radius of the ring. The stretching free energy is 



stretching 



EL 
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in two-dimensions [22, l23[ . The total free energy is then 
given by the sum of the two contributions. 



F = F„ 



F. 



stretching • 



(8) 



Thus, in the inflated phase, the radius of the ring scales 
as 



(9) 



Now, in the inflated phase, {E}) should scale as the 
square of the number of monomers. This implies the 
following scaling 



(10) 



with p = pN . 

In Fig.m we show the variation of {A) /Amax with pres- 
sure p, where A„iax = A^^/16 is the maximum possible 
area. The data points collapse onto one curve when p is 
scaled as p = pN. The data is obtained from exact enu- 
merations of self-avoiding polygons on the square lattice 

M. 



III. 



CONVEX POLYGONS 



In this section we calculate the equilibrium shape and 
area of a convex polygon when p > 0. Convex polygons 
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FIG. 1: The variation of area with pressure p for self avoid- 
ing polygons on a square lattice. Inset: When plotted as a 
function of the scaling variable p = pN, the area curves for 
different values of collapse onto each other. The system 
sizes used are A'' — 50, 60, 70, 80, 90. The data is generated 
from exact enumerations of the polygons on the square lat- 
tice m. 




FIG. 2: A schematic diagram of a convex polygon. Any verti- 
cal or horizontal line (thick dashed lines) intersects the convex 
polygon at either or 2 points. 



are those polygons which have exactly or 2 intersec- 
tions with any vertical or horizontal line drawn through 
the midpoints of the edges of the lattice (see Fig. [5]). We 
calculate the area by determining the shape of the convex 
polygon that minimizes the free energy at fixed perime- 
ter, generalizing the calculation presented in Ref. @. 

The perimeter A'^ of a convex polygon is the same as 
that of its bounding box, which, in general, is a rectan- 
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gle. The equilibrium shape should however be invariant 
about rotations by angle 7r/2. The bounding box of the 
equilibrium shape is thus a square of side N/A. We can 
now calculate the shape in the first quadrant, obtaining 
the shapes in other quadrants by symmetry. 

Consider a coarse grained shape y{x) in the first quad- 
rant with endpoints at (0, iV/8) and (iV/8,0). The free 
energy functional for this curve y{x) can be written as 



N/8 



P_ 

N 



N/8 



dxy, (11) 



where a{y') is the free energy per unit length associated 
with a slope y' and p is the scaled pressure. The shape is 
then obtained obtained from Eq. (fTT|) through the Euler 
Lagrange equation (Wulff construction) [25l] . 



d d 
dx dy' 



P_ 

' N' 



(12) 



The free energy can be calculated using a simple com- 
binatorial argument. Consider all possible paths start- 
ing from {0,y) and ending at (a;,0) with only right- 
ward and downward steps. The weight of a path is 
exp(— J6). where b is the number of bends. When 
x,y ^ 1, the weighted sum of these paths will be equal 
to exp[— + y^cr{y')], where y' = —y/x. 

Let C{x, y) be the sum of weighted walks constructed 
as above. Such walks may be enumerated by splitting 
the path into sequences of rightward and downward steps 
and associating the bending energy term to a sequence 
of downward (y) steps begun and terminated by a step 
to the right. Then 



(13) 



where Wi is the weight associated with the i*'' step, which 
is given by 

W,^[ail~6y^,o) + Sy^,oi (14) 

and a = e~^'^ . The delta function enforces the constraint 
that the steps taken in the y-direction must total y. The 
summation is over all possible numbers of steps taken in 
the y direction at steps 1,2 .. .x. 

Performing a discrete Laplace transform, we obtain 



1 + auj + auj + ■ ■ ■ 

1 - (1 - a)uj' " 
1 -oj 



(15) 



For large x,y, the inverse Laplace transform can be cal- 
culated using the saddle point approximation. This gives 



(16) 
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FIG. 3: Shape of the convex polygon as obtained from 
Eq. pip for polygons with J = 1. The different shapes cor- 
respond to pressure values p = 1.0, 10.0, 50.0, with the outer 
shape corresponding to largest pressure. 



where 

/(w) = y' ln(w) + ln[l - (1 - a)w] - ln(l - w), (17) 
and uj* satisfies 



dio 



0. 



(18) 



The equilibrium shape can now be obtained from 
Eqs. and The shape satisfies the equation 



cC" -h cciC^" ' ' '(1 - a) - cie^ = 1, (19) 

where X and Y are scaled coordinates defined as X = 
xjN and Y = y/N. The constants of integration are 
fixed by imposing the requirement that the shape should 
be symmetric under the interchange of X and Y, and the 
boundary condition that y{x = N/8) = 0. This gives. 



c = — Ci = 



(1 + e-P/8) - V(l + e-P/8)2 -4e-P/8(l-a) 



2(1 -a) 

and the equilibrium shape can be written as 

ceP^ _ c2eP(^+^)(l -a) + ce^^ = 1. 



(20) 



(21) 



The shapes for different values of the scaled pressure are 
shown in Fig. [3] for a convex polygon with J = 1. 

The area of the convex polygon is obtained from the 
equilibrium shape as 



(A) = AN' 



1/8 



YdX, 



(22) 
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where the factor of 4 corresponds to the four quadrants. 
Domg the integration, we obtain 



{A) = 



16 
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pP/8l 



Li2[c(l-a)e'^/VLi2[c(l-a)] 



(23) 



where Li2 is the dilogarithm function. 

The asymptotic behavior for large p may be calculated 
from Eqs. and When p > 1, the constant c 

can be written as, 



Substituting into Eq. ((23|) . we obtain 



(24) 
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+ 0(e-P/«) 



(25) 

When J = 0, the last term on the right hand side of 
Eq. p5|) is zero and the relation reduces to 



iV2 

(A) = — 
^ ' 16 



1 - 



32^2 



3p2 



+ 0{e-^'^), J = 0. 



(26) 




FIG. 4: A schematic diagram of a column-convex polygon. 
Any vertical line intersects the convex polygon at either or 
2 points. 



IV. 



COLUMN-CONVEX POLYGONS 



In this section, we calculate the shape and area of a 
column-convex polygon when p > 0. Column-convex 
polygons are those polygons which have exactly or 2 
intersections with any vertical line drawn through the 
midpoints of the edges of the lattice. There, is however, 
no such restriction in the horizontal direction (see Fig. [5]). 
We calculate the area by determining the shape of the 
column-convex polygon that minimizes the free energy 
for a fixed perimeter. 

The perimeter of a column-convex polygon has no sim- 
ple relation to its bounding box. We thus introduce a 
chemical potential fj, that couples to the perimeter N. 
Consider a shape y{x) with endpoints at {~/3N, 0) and 
(/37V, 0). The free energy functional for this curve y{x) is 
given by 



^[v{x)]= I dxa{y')y/l + y'^ - — dx y 

J-BN J-BN 



(27) 



where, p is the scaled pressure, p = pN. As before, cr{y') 
represents the free energy per unit length associated with 
a slope y' . The Euler-Lagrange equation then gives the 
shape of the curve equilibrium curve y{x). 

The free energy may be calculated as follows. Consider 
all paths starting from (0,0) to {x,y). Let C{x,y) be the 
weighted sum of all paths. Then, we obtain 



(28) 



Vl,V2,---,Vx 



where Wi is the weight associated with the i*^ step and 
equals 



= /il''■l[a(l-^^«,o) + ^,.,o]/i• 



(29) 



Following the steps outlined previously, we convert the 
(5-function in Eq. (j28p into an integral, obtaining 



C{x,y) 



2ii 



(30) 



where 



1 + (1 - 2a)fi^ + ^(a - l)(p>^ + e-*") 
(1 - ^e^'^)(l - ^e-^") 



(31) 

When y ^ 1, Eq. ([30|) may be evaluated by the saddle 
point method. Denoting y/x by y', we obtain 



[isoy' -\n jj,- In f{i^i, a, So)], (32) 



where sq is the saddle point and is given by, 
d 



dso 



ln/(^, a. So) = iy'. 



(33) 



Substituting the expression for (T{y') into the Euler- 
Lagrange equation (Eq. and using Eq. ([55)) we in- 
tegrate once to obtain an equation for y' . This gives 

,_ /i(a-l)(ce-P^/^-c-ieP^/^) 
^ ~ 1 + (1 - 2a)At2 + n{a - l)(ce-P^/^ + c-^eP^/^) 



/zee 



-px/N 



(34) 



The constant of integration c is fixed by the condition 
that the slope of the equihbrium curve is {y' = 0) at 
X = Q. This gives c = 1. Then we can integrate once 
more to obtain the equation of the equilibrium curve as 



Y{X) = 

--+iln 
P P 



(35) 



{1- ^leP^){l- ^le-P^) 



1 + (1 - a)fi^ + fi{a - l)(eP^ + e-P^) 



As before, X and Y are defined as X ^ x/N and 
Y ~ y/N. The constant of integration ci is fixed by 
the boimdary condition y{x = (3N) = 0. This gives, 



ci = In 



(1 - neP'^){l 



1 + (1 - a)yu2 + ^[a - \){ePP + e-Pl^) ' 



(36) 



The parameter (3 that determines the endpoint of the 
curve is still undetermined. It is chosen to be that (3 that 
minimizes the free energy. The Lagrangian Cq for this 
curve Y{X) is given by substituting Eqs. and ([5^ 
into Eq. 



2(3N I dz 
Jo 



In- 



(1 - A<eP^^)(l 



-pz/3 



1 + (1 - q;)^2 ^(-Q, _ i)(e!5^/3 4- e-P^/') 
The parameter (3 satisfies the equation 



dj3 



= 0. 



This gives 
J) 



2fi[l + fi{a~l)] 



(37) 
In fi . 

(38) 

(39) 



2/x[l + /x(a-l)] 

The chemical potential [i is determined by the con- 
straint that total perimeter is N . This is equivalent to 



dC^ 
d\i 



N 



fi then satisfies the equation, 
1 - M.9 



P 
4 



In- 



In g 



(40) 



(41) 



a 



■An 



ia + b){g+l) + Va^-b^g-l) 
V^2— 52 - (a + 5)(c, + 1) _ ^a^-b^g - 1) ' 

where a and b are given by 

a = 1 + (1 - 2a)/x^ (42) 
b = 2^(a- 1). (43) 

This solves the equilibrium macroscopic shape com- 
pletely. The shapes given by Eq. ([55)) are plotted in Fig-E] 
for column-convex polygons with J = 1.0. 
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FIG. 5: Shape of the convex polygon as obtained from 
Eq. (|35p for polygons with J = 1. The different shapes cor- 
respond to pressure values p = 3.0, 10.0, 50.0, with the outer 
curve corresponding to the largest pressure. Both the X and 
Y axes are scaled by (3. 



We now determine the asymptotic behavior of area 
when p ^ 1. This corresponds to the limit /i ^ 0. In 
this limit, g(/i, J) can be expanded as 

gU,J) = --a + aia-l^ii-aia-lfy^-^OU^). (44) 



and Eq. (|^^ reduces to 



(45) 



Substituting the values of g(/i, J) and /i from Eqs. 
and into Eq. ([57]) . we can obtain the Lagrangian in 
the (^, A^) coordinates to be 

^o(/i, iV) = ^ + ^ [Li2(l - a) - ^] +0(e-P/«). (46) 

The Lagrangian in the (j>, N) coordinates can then be 
obtained by a Legendre transformation as 



N 



(47) 



(48) 



when p 3> 1. The area enclosed by the column-convex 
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polygon is 



A 



~2N 



dp 



327r" 



-rrrr Ll2(l — a) 
pi 



(49) 
(50) 



where the factor 2 in Eq. accounts for the lower 

half plane. Interestingly, Eq. ([50]) is identical to Eq. ([^5]) 
which is the asymptotic area expression for convex poly- 
gons. 
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V. SELF-AVOIDING AND 
SELF-INTERSECTING POLYGONS 

In this section, wc study the asymptotic behavior of 
self-avoiding and self-intersecting polygons. An analytic 
calculation along the lines of those presented for con- 
vex and column-convex polygons is not possible for self- 
avoiding polygons. However, we argue as follows: Con- 
vex polygons have no overhangs and the shape has four 
cusps. Introducing overhangs in one direction gives col- 
umn convex polygons, reducing the number of cusps by 
two. Remarkably, the asymptotic behavior of the area in 
the column-convex case [Eq. ([50] coincides with that for 
convex polygons [Eq. (PS)) ] . It is therefore plausible that 
introducing overhangs in both directions does not affect 
the asymptotic behavior of the area, but merely removes 
the remaining two cusps, yielding a smooth shape. We 
therefore conjecture that the asymptotic behavior of the 
area of self-avoiding polygons is given by 



iV2 



327r^ 



64 

-Li2 1 



P>1. (51) 



For self-intersecting polygons in the inflated phase, it 
is expected that the typical shape of the polygon does 
not intersect itself. Therefore, we argue that the area 
of self-intersecting polygons should also have the same 
asymptotic behavior as in Eq. (j51|l . 

These conjectures may be verified numerically. When 
J = 0, the area of self-avoiding polygons may be obtained 



using exact enumeration data on the square lattice 13 1 . 

24|. 



This data is available for lengths up to TV = 90 
When J 7^ 0, there is no exact enumeration data avail- 
able. Wc therefore resort to Monte Carlo simulations. 
The Monte Carlo algorithm consists of a combination of 
global reflection and inversion moves p6j . The system 
size used was N = 800. 

For self-intersecting polygons, the area may be com- 
puted using exact enumeration methods. We briefly de- 
scribe the algorithm for the case J = 0. The generaliza- 
tion to non-zero J can be found in Ref. [23|. Consider a 
random walk starting from the origin and taking steps in 
one of the four possible directions. For each step in the 
positive (negative) x-direction, we assign a weight e"^^ 



FIG. 6: The asymptotic behavior of the area for self-avoiding 
(SAP) and self- intersecting polygons (SIP). The solid lines 
correspond to the theoretical prediction of Eq. (|51[>. The data 
is in good agreement with Eq. (|51|l . 



(e^^), where y is the ordinate of the walker. The weight 
is then e^"^ for a closed walk enclosing an area A. Let 
Tf^{x,y) be the weighted sum of all A^-step walks from 
(0,0) to {x,y). It obeys the recursion relation, 

TN+iix,y) = e-PyTN{x-l,y) + ePyTNix + l,y) 

+TNix,y^l)+TN{x,y + l), (52) 



with the initial condition 



(53) 



Finally, T^r (0,0) gives the partition function of the self- 
intersecting polygons on a lattice. We used exact enu- 
meration data up to A = 150. 

In the case of all exact enumeration data, for each pres- 
sure point, we extrapolated to A — > oo using finite size 
scaling. The results of the numerical analysis is shown 
in Fig. [5) The numerical data agree very well with the 
theoretical prediction. 



VI. 



CONCLUSIONS 



We now summarize the basic results of this paper. We 
studied the asymptotic behavior of area for large pres- 
sures for a class of polygons. For convex and column- 
convex polygons, we calculated the area using the Wulff 
construction. The asymptotic behavior was observed to 
coincide for both classes of polygons. We therefore con- 
jectured that overhangs are not important in the inflated 
regime, and hence that self avoiding polygons should have 
the same asymptotic behavior. This was verified numer- 
ically. We also showed numerically that self intersections 
were irrelevant in the inflated regime. These results con- 
tinue to remain valid in the presence of a non-zero bend- 
ing rigidity. 

Interestingly, the asymptotic behaviour for continuum 
ring polymers differs from that of lattice polygons. In the 
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continuum, the analogous relation for the area of pres- 
surized rings is asymptotically [l^, 



(A) 
A 



p 



p > 1. 



(54) 



This difference between continuum and lattice models is 
physically sensible in the expanded limit, since curvature 
in the lattice case is concentrated in local regions with 



7r/2 bends but is dclocalizcd along the contour in the 
continuum case. 



While our conjectured result for self-avoiding polygons 
is in good agreement with numerical data, it would be of 
interest to have a mathematically rigorous derivation of 
this result. It may be possible to extend the methods of 
Ref. [2l[ to this problem. 



[1] S. Leibler, R. R. P. Singh, and M. E. Fisher, Phys. Rev. [15] 
Lett. 59, 1989 (1987). [16] 

[2] M. E. Fisher, A. J. Guttmann, and S. G. Whittington, [17] 
J. Phys. A 24, 3095 (1991). [18] 

[3] S. V. M. Satyanarayana and A. Baumgaertner, J. Chem. 

Phys. 121, 4255 (2004). [19] 

[4] E. van Faassen, Physica A 255, 251 (1998). 

[5] V. Privman and N. Svrakic, Directed Models of Poly- [20] 
mers, Interfaces, and Clusters: Scaling and Finite-Size 
Properties (Springer- Verlag, 1989). [21] 
[6] R. Rajesh and D. Dhar, Phys. Rev. E 71, 016130 (2005). 
[7] M. Bousquet-Melou, Discrete Math. 154, 1 (1996). [22] 
[8] E. J. J. van Rensburg, The Statistical Mechanics of In- 
teracting Walks, Polygons, Animals and Vesicles (Oxford [23] 
University Press, 2000). [24] 
[9] K. Y. Lin, J. Phys. A. 24, 2411 (1991). 
[10] M. Bousquet-Melou, J. Phys. A. 25, 1925 (1992). 
[11] M. Bousquet-Melou, J. Phys. A. 25, 1935 (1992). [25] 
[12] R. Brak and A. J. Guttmann, J. Phys. A 23, 4581 (1990). [26] 
[13] L Jensen, J. Phys. A 36, 5731 (2003). 
[14] C. Richard, A. J. Guttmann, and L Jensen, J. Phys. A 
34, L495 (2001). 



J. Cardy, J. Phys. A 34, L665 (2001). 

C. Richard, J. Stat. Phys. 108, 459 (2002). 

J. Rudnick and G. Gaspari, Science 252, 422 (1991). 

G. Gaspari, J. Rudnick, and A. Beldjenna, J. Phys. A 

26, 1 (1993). 

E. Haleva and H. Diamant, Eur. Phys. J. E 19, 461 
(2006). 

M. K. Mitra, G. L Menon, and R. Rajesh, preprint 
arXiy: 0708.3318 

T. Prellberg and A. L. Owczarek, Commun. Math. Phys. 
201, 493 (1999). 

A. C. Maggs, S. Leibler, M. E. Fisher, and C. J. Cama- 

cho, Phys. Rev. A 42, 691 (1990). 

M. E. Fisher, J. Chem. Phys. 44, 616 (1966). 

I. Jensen, Number of sap of given 

perimeter and any area, URL 

|http: //www.ms .unimelb . edu. aii/~iwaii/polygons/series/sqsap_pe 

C. Rottman and M. Wortis, Phys. Rep. 103, 59 (1984). 

N. Madras, A. Orlitsky, and L. A. Shepp, J. Stat. Phys. 

58, 159 (1990). 



